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ABSTRACT 

In core-collapse supernovae, v and V are initially subject to significant self- 
interactions induced by weak neutral currents, which may induce strong-coupling 
effects on the flavor evolution (collective transitions). The interpretation of the 
effects is simplified when self-induced collective transitions are decoupled from 
ordinary matter oscillations, as for the matter density profile that we discuss. In 
this case, approximate analytical tools can be used (pendulum analogy, swap of 
energy spectra). For inverted v mass hierarchy, the sequence of effects involves: 
synchronization, bipolar oscillations, and spectral split. Our simulations shows 
that the main features of these regimes are not altered when passing from sim- 
plified (angle-averaged) treatments to full, multi-angle numerical experiments. 

1. Prologue 

Densely packed inviduals often show a surprising, collective behavior (Fig. [1]). 
Recent developments suggest that neutrinos make no exception, despite the weakness 
of their self-interactions. Effects of u-u forward scattering (via neutral currents) 
may be as important as the known effects of z/g-e" forward scattering in matter (via 
charged currents), provided that the u number density is very high. The dense core of 
exploding Supernovae might provide a possible environment where the flavor evolution 
of z/'s and F's can show, indeed, highly nonlinear and strongly coupled effects. 





Figure 1: In a school of fish, individuals often show a collective behavior. Very dense neutrino gases, 
like those emerging from a core-collapse supernova, might show analogous features in flavor space. 
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Figure 2: Radial profiles of the neutrino self-interaction parameter /i(r) = ^/2Gf (N + N) and of the 
matter-interaction parameter A(r) = \/2Gf N^- adopted in this work, in the range r G [10, 200] km. 



2. Reference supernova model 

Supernova u oscillations are a very important tool to study astrophysical pro- 
cesses and to better understand u properties After leaving the neutrinosphere, 
u and u undergo flavor oscillations triggered by vacuum mass-mixing parameters 
and by ordinary (MSW) matter effects. Besides, in the flrst few hundred kilometers 
neutrino-neutrino interactions may induce additional important effects (depending on 
the neutrino mass hierarchy). Self-interaction effects are expected to be non negli- 
gible when yu(r) ~ cu, where fi = ^/2GF{N^{r) + N,y{r)) is the potential associated 
to the u + T' background [analogous to the MSW potential A = \/2GFNe'{r)], while 
u = Am'^/E is the largest vacuum oscillation frequency. We neglect the smallest 
mass squares difference Sm"^ = m\ — m\ <^ Aim? = [mg — m^2l) consider a 2i/ 
mixing scenario governed by Am^ and the mixing angle Oi^ {Am? = 10~^ eV^ and 
sin^^is = 10~^ for reference). In the supernova context, and Ur (shortly, u^) be- 
have similarly, and we can generically consider two-neutrino z/g oscillations as a 
reasonable approximation. 

Figure 1 shows the radial proflles of the matter potential A(r) and of the neutrino 
potential /u(r), and the approximate ranges where different collective effects occur: 
synchronization, bipolar oscillation and spectral split. The nonlinearity of the self 
interactions induce collective transitions for small r, well before the ordinary MSW 
resonance, allowing a clear interpretation of the numerical simulations. However, for 
matter proflles different from ours (shallow electron density proflles®), the MSW 
effects can be already operative around O(IOO) km, in which case (not typical, and 
not considered here) they are entangled to the collective ones in a complicated way. 



We adopt normalized thermal spectra with [Ee) = 10 MeV, (Ee) = 15 MeV, and 
(Ere) = (Ex) = 24 MeV for Ug, Vg, and Vx, respectively. The emission geometry 
is based on the so called "bulb model" ^ with spherical symmetry: neutrinos are 
assumed to be half-isotropically emitted from the neutrinosphere. Along any radial 
trajectory there is, therefore, a cylindrical symmetry. As cylindrical variables one 
can choose the distance form the supernova center r, and the angle d between two 
interacting neutrino trajectories. If the dependence on d is integrated out, one speaks 
of "single-angle" approximation, while the general situation of variable d is dubbed 
"multi-angle" case. The numerical simulation in the multi-angle case is extremely 
challenging, since it typically requires the solution of a large system (~10*^, M > 5) 
of coupled non-linear equations, after discretization of the cylindrical coordinates. 

3. Equations of motion, pendulum analogy, and spectral split 

The propagation of neutrinos of given energy E is studied through the Liouville 
equation for the 2x2 neutrino density matrix in flavor basis. By expanding it on 
the Pauli and the identity matrices, the equations of motion can be expressed in 
terms of two flavor polarization vectors, P(-E') for any neutrino and P(-E) for any 
antineutrino. By introducing a vector B that depends on the mixing angle ^13, and 
a vector D = J — J that is the difference between the integral over the energy of P 
and P, the equations of motion can be written as, 

P = (+wB + Az + /iD) X P , (1) 
F = (-cuB + Az + /iD) X P , (2) 

see^^ and references therein. In the general case, the polarization vectors depend also 
on the neutrino emission angle 9q at the neutrinosphere (the neutrino intersection 
angle d can be expressed in terms of r and of ^o)- The survival probability Pee is 
a function of the flavor polarization z-component. Pee = 1/2(1 + / PI)., where the i 
and / refer to the initial and final state, respectively (analogously for Vf,). 

Collective effects show up by aligning such polarization vectors (in flavor space) 
close to each other. In the alignment approximation, the equations of motion for P(-E') 
and P{E) can be reduced to collective ones describing a classical, gyroscopic pendu- 
lum: a spherical pendulum of unit length in a constant gravity field, characterized by 
a point-like ma ssive bob spinning around the pendulum axis with constant angular 
momentum . The pendulum inertia is inversely proportional to /i(r), while its 
angular momentum depends on the difference of the integrated polarization vectors 
J and J, see®. The motion of a spherical pendulum is, in general, a combination of 
a precession and a nutation. 

In the case of normal neutrino mass hierarchy, the pendulum starts close to the 
stable, downward position and stays close to it, as yU slowly decreases collective effect 
gradually vanish. In the inverted hierarchy case, the pendulum starts close to the 




Figure 3: Single- angle simulation in inverted hi- 
erarchy: Pz (neutrinos) and Pz (antineutrinos) 
as a function of radius, for five energy values. 



Figure 4: Multi-angle simulation in inverted hi- 
erarchy: Pz (neutrinos) and Pz (antineutrinos) 
as a function of radius, for five energy values. 



"unstable," upward position, being slightly tilted by an angle of 0(^13). At small r, 
when n is large (small pendulum inertia), the bob spin dominates and the pendulum 
remains precessingin the upward position as a "sleeping top" a situation named 
synchronization As /i decreases with r, the pendulum inertia increases and, 

unavoidably, for any ^^13 7^ 0, the pendulum fall occurs with subsequent nutations, 
the so called bipolar oscillations. The increase of the pendulum inertia with r reduces 
the amplitude of the nutations, and bipolar oscillations are expected to vanish when 
self-interaction and vacuum effects are of the same size. 

While the bipolar regime comes to an end, self-interaction effects do not completely 
vanish, and a spectral split builds up: a "stepwise swap" between the and energy 
spectra. The neutrino swapping can be explained by the conservation of the pendulum 
"energy" and of the lepton number ^ . The lepton number conservation is related to 
the constancy of = Jz~ J z-, that is a direct consequence of the equation of motion. 
For a more detailed description of the pendulum analogy in our reference model, the 
reader is referred to our work ^ and references therein. 



4. Single- and multi-angle simulations: stability of results 

Figures 3 and 4 show the z-component of P and P, as a function of r at different 
E values, for single- and multi-angle simulations, respectively. Bipolar oscillations 
start after a synchronization plateau, with equal periods for both v and V at any 
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Figure 5: Single- angle simulation in inverted hi- 
erarchy: modulus and z-component of J and J. 



Figure 6: Multi-angle simulation in inverted hi- 
erarchy: modulus and z-component of J and J. 



energy, confirming the appearance of collective features. Indeed, the behavior of each 
Pz and Pz depends essentially on its energy. For neutrinos. Figure 3, the spectral split 
(inversion of Pz) starts at a critical energy Ec — 7 MeV: the curve relative to E < Ec 
ends up at the same initial value (Pee = 1), while the curves for E > Ec show the Pz 
inversion (Pee = 0). Neutrinos with an energy of ~ 19 MeV do not oscillate much, 
because this is roughly the energy for which the initial and fluxes are equal in 
our scenario. For 77, all curves show almost complete polarization reversal, except at 
very small energies (of few MeV, not shown). Multi-angle simulations (Fig. 4) are 
similar, although with somewhat damped bipolar oscillations. 

Figures [5] and [6] show the evolution of the global polarization vectors (modulus J 
and 2-component Jz) for neutrinos and antineutrinos, in the single- and multi-angle 
cases. The behavior of these vectors can be related to the gyroscopic pendulum mo- 
tion. At the beginning, in the synchronized regime, all the polarization vectors are 
aligned so that J = Jz and J = Jz'- the pendulum just spins in the upward position 
without falling. Around ~ 70 km the pendulum falls for the first time and nuta- 
tions appear. The nutation amplitude gradually decreases and bipolar oscillations 
eventually vanish for r ~ 100 km. At the same time, the spectral split builds up: 
antineutrinos tend to completely reverse their polarization, while this happens only 
partially for neutrinos. As said before, also for antineutrinos there is a partial swap 
of the spectra for ~ 4 MeV. From Figure [6] it appears that bipolar oscillations of 
J and J are largely smeared out in the multi-angle case. The bipolar regime starts 
somewhat later with respect to the single-angle case, since neutrino-neutrino inter- 
action angles can be larger than the (single-angle) average one, leading to stronger 
self-interaction effects, that force the system in synchronized mode slightly longer. 
However, just as in the single-angle case, the spectral split builds up, Jz gets finally 
reversed, while the difference Dz = Jz — Jz remains constant. 



Final fluxes in inverted hierarcliy (single-angle) 
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Figure 7: Single-angle simulation in inverted hierarcliy: final fluxes (at r = 200 km, in arbitrary 
units) for different neutrino species as a function of energy. Initial fluxes are shown as dotted lines. 

Figures [7] and [S] show the final v and V fiuxes, in the single- and multi-angle 
simulations. The v clearly show the spectral split effect and the corresponding sudden 
swap of Vf. and fiuxes above E^. — 1 MeV. In the right panel of Figure [71 the 
final V spectra are basically completely swapped with respect to the initial ones, 
except at very low energies, where there appears a minor V spectral split. This 
phenomenon can be related to the loss of J and of \Jz\ Also in the multi-angle 
case of Figure [8], the v spectral swap at E > Ec — 7 MeV is rather evident, although 
less sharp with respect to the single-angle case, while the minor feature associated to 
the "antineutrino spectral split" is largely smeared out. We conclude that at least 
the I' spectral split (left panel of either Fig. 7 or 8) provides a robust and potentially 
observable collective feature emerging in numerical experiments for inverted hierarchy. 



Final fluxes in inverted hierarchy (multi-angle) 
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Figure 8: As in Fig. 7, but for multi-angle simulations. 



5. Summary 



We have studied supernova neutrino oscillations in a model where the collective 
flavor transitions (synchronization, bipolar oscillations, and spectral split) are well 
separated from later, ordinary MSW effects. We have performed numerical simu- 
lations in both single- and multi-angle cases, using continuous energy spectra with 
significant v-V and i^e-^x asymmetry. The results of the single- angle simulation can 
be largely understood by means of an analogy with a classical gyroscopic pendulum. 
The main observable effect appears to be the swap of final-state energy spectra, for 
inverted hierarchy, at a critical energy dictated by lepton number conservation. In 
the multi-angle simulation, details of self-interaction effects can change, but the spec- 
tral split remains a robust, observable feature. In this sense, averaging over neutrino 
trajectories does not alter the main effect of the self interactions. From the point 
of view of neutrino parameters, collective flavor oscillations in supernovae could be 
instrumental in identifying the inverse neutrino mass hierarchy, even for tiny ^13. ^ 

6. Epilogue 

Supernova u and V flavor polarization vectors can perform elaborate and collective 
"dances" (precession, nutations) in flavor space, at least for the first O(IOO) km, in 
the case of inverted mass hierarchy. Many aspects of this behavior, however, seem to 
be precluded to experimental observations, with the possible exception of a robust, 
finale-state feature: the v spectral split. According to a calculable energy threshold, 
supernova v might then proceed to the Earth with their original flavor (say, z/g), or 
with the complementary one (say, z/^,), just as an initially coherent school of fish 
(Fig. [1]) may finally branch out in some circumstances (Fig. [9]). 




Figure 9: A school of fish branching out in two different directions. Analogously, supernova neutrino 
polarization vectors might split up in flavor space, due to self-interaction effects. 
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